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(54) METHOD AND SYSTEM FOR OPTIMIZING PLANT UTILITY 



(57) Optimum operation parameters of a plant are 
sought, the plant including a plurality of first plant ele- 
ments B1 to B4 having linear input/ouput characteristic 
that receive first physical quantities (f1 to f4) and output 
second physical quantities (s1 to s4), and a plurality of 
second plant elements (T1 to T5) that receive the sec- 
ond physical quantities (s1 to s4), i.e., the outputs from 
the first plant elements (B1 to B4), and output third 
physical quantities, the second plant elements including 
non-salient characteristic portions. For this purpose, the 



group of second physical quantities and the group of 
third physical quantities are processed twice by utilizing 
a dynamic planning method, and are processed one 
time by utilizing a linear planning method, in order to 
find a solution for optimizing the output of the group of 
second plant elements for every operation condition. 
The operation condition in which the total cost is mini- 
mized for the first physical quantities is determined to be 
the optimum solution. 
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Description 

[Technical Field] 

5 The present invention relates to a plant utility optimizing method for optimizing an operation state of a plant system 

while satisfying various demands for the plant system and to an optimizing system. 

[Background Art] 

10 In a factory having, for example, normally used steam turbine generators, an optimal operation system for a power 
generation plant operates the power generation plant (hereinafter referred to as a BTG plant) so as to supply energy 
which is required in a manufacturing process to power and steam loads at a minimal cost. 

FIG. 18 shows a conceptual diagram of the BTG plant and FIG. 19 shows a practical system arrangement of the 
BTG plant. 

15 As shown in FIG. 1 8, a plurality (Xm) of steam turbines T (hereinafter referred to simply as the turbines) are oper- 
ated by the steam generated by steam generation boilers B (hereinafter referred to as the boilers) using fuel F such as 
the oil, gas, coal, recovered black liquor, etc. The respective turbine T drives a corresponding generator G to generate 
electric power. The electric power is supplied to electric power loads in the associated plant The steam bled from the 
turbine T is supplied to the steam loads in the plant. In this case, rf only an electric power generated from the generators 

20 is not sufficient to meet the electric load demand in the plant, a commercial power is supplied to the loads. 

The system arrangement of the BTG plant comprises, as shown in FIG. 19 for instance, boilers B1, B2 using the 
well-controlled fuel such as the oil, gas, etc., boilers B3, B4 using the fuel of better economic efficiency such as the coal, 
black liquor, etc., a high-pressure steam header SL supplied with the steam generated from the boilers B1 , B2, a high- 
pressure steam header SH supplied with the steam generated from the boilers B3, B4, steam turbines T1 , T2 operated 

25 by the steam pressure coming from the high-pressure vapor header SH and driving the generator G, steam turbines T3, 
T4 operated by the steam pressure coming from the high-pressure steam header SL and driving the generator G, a 
supply header DH for supplying the steam which is bled from the steam turbines T1 , T2 to the associated loads, and a 
supply header DLfor supplying the steam which is bled from the steam turbines T3. T4 to the associated loads. Here, 
the pressure of the high-pressure steam header SH is higher than that of the high-pressure steam header SL, while the 

30 pressure of the steam supply header DH is higher than that of the steam supply header DL. 

The high-pressure steam and low-pressure steam supplied to the headers DH and DL are supplied, respectively, 
to the steam loads H and L in the plant. 

A pressure-reducing valve RPV1 is connected between the high-pressure steam header SH and the high-pressure 
steam header SL and a pressure-reducing valve is connected between the headers DH and DL. 

35 in the BTG plant, the fuel amount/steam generation amount characteristic of the boiler shows a gentle, near-flat 
curve as shown in FIG. 20. According to the boiler characteristic as shown in FIG. 20, the evaporation amount from the 
boiler is found to be substantially in proportion to the fuel amount in an upper/lower limit operation range. 

On the other hand, the steam amount/generation power output characteristic of the turbine has a non-convex char- 
acteristic called a valve loop characteristic (valve point characteristic), the characteristic curve having several recesses. 

40 The non-convex characteristic is due to a fall in a shaft output resulting from a pressure drop produced when the steam 
input valve is half-opened. In the characteristic curve proposed by the apparatus manufacturer as the characteristic 
curve of the turbine, a plurality of peak points are normally connected together to provide an envelope, an aspect which 
is different from a real aspect shown in FIG. 21 . That is, according to the turbine characteristic as shown in FIG. 21 it 
has been found that the generation output power differs in the upper/lower limit range in terms of the valve points and 

45 their intervening recesses. 

Incidentally, at the optimized operation of the BTG plant the steam distribution of the boiler and turbine and that of 
the commercial power are varied so that the operation is performed at a minimal cost. This effect will be explained with 
reference to FIGS. 22 and 23. 

Let it be assumed that, as shown in FIG. 22, two turbines of the same characteristic are operated and, at this time, 
so the steam load is at 1 00 T/H. 

Here, the load steam, being equally distributed to the two turbines, becomes a 50 T/H each and, at this time, the 
generated power amount becomes 

200 x 50 + 200 x 50 = 20000 KW 

55 

as shown in FIG. 23. 

If, on the other hand, the load steam is distributed to the two turbines, one at 40 T/H and the other at 60 T/H, the 
generated power output becomes 210 x 40 + 218 x 60 = 21480 KW. From this it has been found that 7.4 % more output 
can be generated by varying the distribution of the steam to the turbines. 
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In the optimizing operation system of the BTG plant, therefore, the BTG plant is modeled with a mathematical for- 
mula and a cost minimum is found from boiler/turbine distribution combinations which can satisfy the electric 
power/steam loads. 

Given that a target total energy cost is expressed with a function, the BTG plant model with a mathematically for- 
s matted distribution amount at which the function becomes minimal (or maximal) is found, as a mathematical solution, 
with a high-speed computer. By operating the BTG plant to reach this distribution amount it is possible to operate, at a 
minimal cost at all times, the BTG plant which would be complex and impossible to the human. 

FIG. 26 is a block diagram showing one form of an optimal operation system of this kind. In FIG. 26, the optimal 
operation system of the BTG plant comprises a plant model 32 and a total energy cost minimizing point calculating see- 
to tion31. 

The characteristics and operation restrictions of the turbines, boilers, pressure-reducing valves, etc., as the constit- 
uent elements of the BTG plant, as well as the electric power balance and steam balance at each system part, are pro- 
vided from the plant model 32 to the total energy cost minimizing point calculating section 31 . 

By giving the actual load amount (electric power and steam) as the plant state of the plant model 32 to the total 
15 energy cost minimizing point calculating section 31 an optimal distribution is found with the use of a linear programming 
(LP) or a nonlinear programming (NLP). 

The linear or nonlinear programming conventionally used thus far has advantages and disadvantages. 

That is, the mathematically formulated model as a target in the linear programming can solve only all linear (linear 
equation) problem and it is necessary, therefore, to apply the plant model, being nonlinear in actual practice, through a 
20 linear approximation. For this reason, it is not possible to avoid a lowering in the accuracy of a solution through such 
approximation, so that the cost lowering effect is small. 

The non-linear programming can solve the multi -dimensional program and, therefore, handle the nonlinear plant 
model but, in the case where the problem to be handled is non-linear, it is possible to obtain only a local optimization 
solution with the attendant drawback. 
25 In the case where, for example, a valve point characteristic of the turbine is non-convex and multi-peaked, it is pos- 
sible to. for a local optimization point and global optimization point as shown as the conceptual diagram, obtain only an 
optimal point (local optimization solution) in the vicinity of a calculated initial value so that no true optimal point (global 
optimal solution) can be obtained. It is because the object function is multi-peaked. 

Further, in order to obtain a solution as near to a true optimal point as possible with the use of the nonlinear pro- 
30 gramming, the initial value and search method have to be modified in accordance with the problem. However, the 
search logic beoomes complex when aiming at an optimal solution and lots of time and labor are taken in tuning calcu- 
lation parameters. 

For the search of the optimal point with the use of the nonlinear programming (NLP) it is necessary to usually deter- 
mine the initial value and search range. 
35 The initial value is used to determine from where the calculation is started. The method of how to determine the 
initial value is often by, for example, 

i) using the current operation state, 

ii) getting the respective characteristic functions all linearly approximate and utilizing a solution obtained by the lin- 
40 ear programming (LP). 

Further, for the Lagrande multiplier method-plus-conjugate gradient method often used as the NLP alogrithm, the 
search range has to be tuned to the problem so that a single peak (trough) is obtained in the search range. 

FIG. 25 shows a relation among the linear, non-linear and true optimal resolutions. 
45 In FIG. 25, the optimal solution by the LP is located at a T Dj point. 

Further, when the initial value is started from a \ Kj point by the NLP, it is possible to obtain, as an optimal solution, 
a minimal point I" AJ in the vicinity of the search done. This, being greater than a true optimal point \ TJ, provides a local 
optimal point. 

A point T Bj is one shifted by L2 from the local optimal point I" AJ at the search of a better point. In the case where 
so the T BJ point is located in a position lower than the I" AJ point, the optimal point is shifted to the I" Bj point. And a local 
optimal point [ TJ in the vicinity of the search is obtained as an optimal solution as in the case of the I" AJ point. 

However, a point I" CJ spaced apart by L1 is located in a position higher than the I" AJ point, so that there is no shift 
of the optimal point. 

In the search of the optimal point with the use of the NLP method there is no right determination of the search range 
55 L and hence no guarantee to reach a true optimal point [ TJ. 

Since the non-linear programming is based on the assumption that the problem to be dealt with is on the convex 
characteristic, it cannot be used unless it is known whether there is one peak (single peak) and wether the position of 
the peak is located. 
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It is found that no local optimal point is searched on a multi-peak problem containing a plurality of non-convex char- 
acteristics as in the case of the valve characteristics of the turbine. 

In the conventional optimal system, therefore, the operation is done at the local optimal point and it may be said 
that there is room for efficiency improvement up to the local optimal point. 

s For the nonlinear programming method it has been difficult to tune calculation parameters such as the initial value 
and a way of how to determine the search range and it has been very difficult, for the user with no expert's knowledge 
on the optimization computation, to do the system changes of the target energy plant and to do the system mainte- 
nance such as the correction of a plant model corresponding to apparatus changes 

,„ *!! f ' ndin9 ° f *" 0pt ' mal distribution at * e tot al energy cost optimization point calculating section it has been pos- 

w sible for the conventional plant utility optimization system to obtain a linear solution of poor accuracy only or a nonlinear 
solution allowing a local optimization solution only and it has been difficult to. for the nonlinear solution, tune the calcu- 
lation parameters such as the initial value and search range. 

A first object of the present invention is to provide a plant utility optimization method which can satisfy various 
demands for a plant system and find an optimal operation state reflecting non-convex characteristics possessed by 

15 plant elements. 3 

A second object of the present invention is to provide a plant optimization system which can ensure the obtainment 
of a solution of higher accuracy than in a conventional system and can obviate the need to tune calculation parameters 
involved and. in addition, enables the user with no expert s knowledge on optimization calculation to do the system 
20 to 2plratus l chL 9 ngeT er9y **** * * SyStGm maintenance such as the <w«tton °* a P'ant model corresponding 

[Disclosure of Invention] 

A utility optimizing method of the present invention is provided for searching optimal operation parameters of a 
25 plant including a plurality of first plant elements having a linear input/output characteristic with first physical amounts as 

3 h a T 3 1 Se S??„ d PhySi °!! am ° UntS 38 an OUtpu1 a p,urality d second P |ant element s ^ving an input/output 
characteristic with the second physical amounts as an input, conesponding to the output of the first plant element and 
third physical amounts as an output, the input/output characteristic containing a non-convex characteristic portion the 
method comprising: 



a first processing step of finding a distribution of the second physical amounts to the second plant elements by a 
dynamic programming to maximize the sum total of the third physical amounts corresponding to the sum total of 
the second physical amounts in a predetermined range satisfying a first demand of the second physical amounts 
in a first group and output from the second plant elements associated with the first group- 

35 f second processing step of finding a distribution of the second physical amounts to the second plant elements by 
thedynamic programming on the basis of a set of the sum total of the second physical amounts, sum total of the 
hird physical amounts and distribution of the second physical amounts to the second plant elements obtained by 
the first processing step to maximize the sum total of the third physical amounts corresponding to the sum total of 
the second physical amounts to the second plant elements and output from the respective second plant elements- 

40 and ' 

a third processing step of finding by a linear programming or nonlinear programming, the first physical amounts 
input to the second plant elements and evaluation values, necessary for supplying the first physical amounts with 
r^pect to each set of the sum total of the second physical amounts, sum total of the third physical amounts and 
distribution of the second physical amounts to the second plant elements obtained at the second processing step. 

Further, a utility optimizing method of the present invention is provided for searching optimal operation parameters 
of a plant including a plurality of first plant elements having a linear input/output characteristic with first physical 
amounts as an input and second physical amounts as an output and a plurality of second plant elements having a linear 
input/output characteristic with the second physical amounts as an input, corresponding to the output of the first plant 
elements, and third physical amounts as an output, the input/output characteristic containing a non-convex character- 
istic portion, the method comprising: 

a first processing step of finding a distribution of the second physical amounts to the second plant elements by a 
dynamic programming to maximize the sum total of the third physical amounts corresponding to the sum total of 
the second physical amounts in a predetermined range satisfying a first demand for a group of such second phys- 
ical amounts and output from the second plant elements associated with the group- 

a second processing step of finding a distribution of the second physical amounts to the second plant elements by 
a dynamic programming on the basis of a set of the sum total of the second physical amounts, sum total of the thiid 
physical amounts and distribution of the second physical amounts to the second plant elements obtained by the 
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first processing step at each second physical amount in a group to maximize the sum total of the third physical 
amounts corresponding to the sum total of the second physical amounts to all the second plant elements and out- 
put from said the second plant elements; 

a third processing step of finding, by a linear or a nonlinear programming, the first physical amounts input to the 
5 second plant elements and first evaluation values, necessary to supply the first physical amounts, with respect to 
the respective set of the sum total of the second physical amounts, sum total of the third physical amounts and dis- 
tribution of the second physical amounts to the second plant elements obtained by the second processing step; 
a fourth processing step of. on the basis of a respective set of the sum total of the third physical amounts and first 
evaluation values, necessary for the first physical amounts, calculating second evaluation values required to 
to receive a power which corresponds to a shortage in a predetermined second demand for the sum total of the third 
physical amounts to find a total evaluation value for each set and deciding, as an optimal solution, each value 
belonging to a set in which the total evaluation value is minimized. 

Further, a utility optimizing method of the present invention is provided for searching optimal operation parameters 
15 of a power generation plant including a plurality of boilers, a plurality of turbines operated by steams generated from the 
boilers, and generators driven by the respective turbines to generate electric power in which the electric power is sup- 
plied to electric power loads while the steams bled from the respective turbines are supplied to steam loads via a plu- 
rality of turbine discharge systems, the method comprising: 

20 a first processing step of finding a steam flow amount distribution to the respective turbines by a dynamic program- 
ming to maximizing the sum total of generation output powers of the turbines corresponding to respective steam 
flow amounts in a predetermined range satisfying a steam demand to the respective turbine discharge system and 
connected to the same turbine discharge system; 

a second processing step of, on the basis of a set of the steam flow amount, generation output power and steam 
25 flow amount distribution obtained by the first processing step at each turbine discharge system, finding a steam 
flow amount distribution to the respective turbines by a dynamic programming to maximize a total generation output 
power of all the turbines corresponding to the total input steam flow amounts to all the turbines; and 
a third processing step of finding a supply fuel distribution to the respective boilers and total fuel cost by a linear or 
a nonlinear programming with respect to each set of the total input steam flow amount, total generation output 
30 power and steam flow amount distribution to the respective turbines. 

Further, a utility optimizing method of the present invention is provided for searching optimal operation parameters 
of a power generation plant including a plurality of boilers, a plurality of turbines operated by steams generated from the 
boilers, and generators driven by the respective turbines to generate electric power in which the electric power is sup- 
35 plied to electric power loads while the steams bled from the respective turbines are supplied to steam loads via a plu- 
rality of turbine discharge systems, comprising the steps of 

a first processing step of finding a steam flow amount distribution to the respective turbines by a dynamic program- 
ming to maximize the sum total of generation output powers of the turbines corresponding to the steam flow 
40 amounts in a predetermined range satisfying a steam demand to the respective turbine discharge system and con- 
nected to the same turbine discharge system; 

a second processing step of finding a steam flow amount distribution to the respective turbines by the dynamic pro- 
gramming on the basis of a set of the steam flow amount, generation output power and steam flow amount distri- 
bution obtained by the first processing step at each turbine discharge system to maximize a total generation output 
45 power of the turbines corresponding to a total input steam flow amount to all the turbines; 

a third processing step of finding, by a linear or nonlinear programming, a supply fuel distribution to the respective 
boilers and total fuel cost with respect to each set of the total input steam flow amount, total generation output 
power and steam flow amount distribution to the respective turbines; and 

a fourth processing step of, on the basis of each set of a total generation output power and total fuel cost obtained 
so by the third processing step, calculating an external cost required to receive a power which corresponds to a short- 
age in a predetermined power demand for the total generation output power from an outside source to find a total 
cost for each set and deciding, as an optimal solution, each value belonging to a set in which the total cost is min- 
imized. 

55 Further, in a power generation plant including a plurality of boilers, a plurality of turbines operated by steams gen- 
erated from the boilers and generators driven by the respective turbines to generate electric power wherein the electric 
power is supplied to electric power loads while the steams bled from the respective turbines are supplied to the steam 
loads, a plant utility optimizing system of the present invention comprising: 
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a plant model prepared by separating into a turbine-side and a boiler-side model; 

a turbine output maximum point calculating section for deciding such a steam distribution of the turbines as to max- 
imizing a total of outputs of the respective turbines for each supply steam amount given from the turbine-side model 
in the plant model; 

5 a fuel cost minimizing point calculating section for deciding such a steam distribution amount of each boiler as to 

generate at a minimal cost, the sum total of per pressure supply steam amounts given from the boiler-side model 
at each point found by the turbine output maximum point calculating section; and 

a total energy cost minimizing point calculating section for searching a total energy cost minimum operation point 
to electric loads and steam loads on the basis of a steam distrfoution amount decided by the fuel cost minimizing 

10 point calculating section. Further, in a power generation plant including a plurality of boilers, a plurality of turbines 
operated by steams generated from the boilers and generators driven by the respective turbines to generate elec- 
tric power in which the electric power is supplied to electric power loads while steam bled from the respective tur- 
bine is supplied to steam loads, a plant utility optimizing system of the present invention comprising: 
a plant model prepared by separating into a turbine-side and a boiler-side model; 

is a turbine output maximum point calculating section for deciding such a steam distribution of the turbines as to max- 
imize a total of outputs of the respective turbines for each supply steam amount given from the turbine-side model; 
a fuel cost minimizing point calculating section for deciding such a steam distribution amount of the respective boil- 
ers as to generate, at a minimal cost, the sum total of per pressure supply steam amounts given from the boiler- 
side model at each point found by the turbine output maximum point calculating section; 

20 a total energy cost minimizing point calculating section for searching a total energy cost minimum operation point 
to electric power loads and steam loads on the basis of a steam distribution amount of the respective boilers 
decided by the fuel cost minimizing point calculating section; 

a plant model preparing section for preparing/editing system graphics for modifying or correcting plant models used 
under the calculation of the turbine output maximum point calculating section and fuel cost minimizing point calcu- 
25 latmg section and system models described with a mathematical expression; and 

a data generator for generating data necessary for optimal calculation from the model prepared by the plant model 
preparing section and for giving the data to the turbine-and boiler-side models in the plant model. 

According to the present invention, in the searching of optimal operation parameters in a plant (for example, the tur- 
bine-type power generation system) including a plurality of first plant elements having a linear input/output characteris- 
tic with first physical amounts (for example, fuel) as an input and second physical amounts (for example, steam 
amounts) as an output and a plurality of second plant elements having an input/output characteristic with the second 
physical amounts as an input, corresponding to the output of the first plant elements, and the third physical amounts 
(for example, output power) as an output, the latter input/output characteristic containing a non-convex characteristic 
portion, processing is done twice through the utilization of the dynamic programming allowing the non-convex charac- 
teristic to be reflected. Once the first demand (for example, steam demand), that is, a demand for the second physical 
amounts, has been satisfied, sets of the sum total of the second physical amounts, distribution amount of this sum total 
to all the second plant elements and sum total of the third physical amounts output from all the second plant elements 
are found, noting that the total amount of various second physical amounts is given here. 

For each set, the first physical amounts input to the respective second elements and first cost required to supply 
these first physical amounts are found by a linear or a nonlinear programming on the basis of the values in the found set. 

Finally, in the case where, for each set, the sum total of the third physical amounts is smaller than a second demand 
(for example, an energy demand), a second cost (for example, external cost) required to receive its shortage from an 
outside is calculated on the basis of the values in the set of the sum total of the third physical amounts and first cost (for 
example, a total fuel cost) required for the first physical amounts thus obtained to find a total cost for each set and a 
respective value (for example, the sum total of the second physical amounts, distribution of the second physical 
amounts to the second plant elements, sum total of the third physical amounts, first physical amounts input to the 
respective second plant elements, first cost required to supply the first physical amounts, total cost, etc.,) belonging to 
a set, that is, a set in which that total cost is minimized, is decided as an optimal solution. 

According to the plant utility optimizing system of the present invention, it is possible to decide such a steam distri- 
bution (point) of the respective turbines as to maximize a total of the outputs of the respective turbines for each supply 
steam amount given from the turbine-side model in the plant model and it is also possible to decide such a steam dis- 
tribution amount of the respective boilers as to generate, at a minimal cost, the sum total of the per pressure supply 
steam amounts given from the boiler-side model for each point. By searching a total energy cost minimum operation 
point (including a power buying cost) for the electric loads and above-mentioned loads it is possible to obtain a true opti- 
mal point and largely lower a cost involved. 

In addition to the above-mentioned functions, it is also possible to, on the plant model used for calculating the tur- 
bine output maximum point and fuel cost minimizing point, prepare/edit the system graphics and system model 
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described with the mathematical expression and generate, from this model, data necessary for optimal calculation. By 
doing so, it is possible to readily modify or correct the model 

[Brief Description of Drawings] 

5 

FIG. 1 diagrammatically shows a plant energy system according to a first embodiment of the present invention; 
FIG. 2 is a view showing one example showing a schematic flow relating to a plant utility optimizing system accord- 
ing to the first embodiment of the present invention; 

FIG. 3 is a view showing a major section of the plant system according to the first embodiment; 
w FIG. 4 is a schematic view showing a major section of the plant utility optimizing system according to the first 
embodiment; 

FIG. 5 is a flow chart showing one example of the operation of a fluid demand system optimizing section; 
FIG. 6 is a flow chart showing the example of the operation of the fluid demand system optimizing section; 
FIG. 7 is a view showing a relation of an optimal steam distribution and output power to each steam amount of an 
15 L3 group; 

FIG. 8 is a view showing a relation of an optimal steam distribution and output power to each steam amount of an 
L4 group; 

FIG. 9 is a view showing a relation of an optimal steam distribution and output power to each steam amount of an 
L3 group; 

20 FIG. 1 0 is a view showing a relation of an optimal steam distribution and output power to each steam amount of an 
L4 group; 

FIG. 1 1 is a view showing a relation of an optimal steam distribution and output power to each steam amount of an 
L3 group; 

FIG. 12 is a view showing a relation of an optimal steam distribution and output power to an L4 group; 
25 FIG. 1 3 is a flow chart showing one example of the operation of a fluid input system optimizing section; 

FIG. 14 is a view showing a relation of an optimal steam distribution and total output power to a total input steam 
amount; 

FIG. 15 is a view showing a relation of an optimal steam distribution, total output power, fuel distribution and total 
fuel cost to a total input steam amount; 
30 FIG. 16 is a view showing respective optimal solution candidates; 

FIG. 1 7 is a block diagram showing an optimal operation system of a power generation plant according to a second 

embodiment of the present invention; 

FIG. 18 is a conceptual view showing a BTG plant; 

FIG. 19 is a view showing an example of a system arrangement of a BTG plant; 
35 FIG. 20 is a characteristic curve showing a relation of the fuel of boilers used in the BTG plant to a steam generation 
amount; 

FIG. 21 is a characteristic curve showing a relation of a steam generation amount of the turbines used in the BTG 
plant to a generation power output; 

FIG. 22 is a view showing an example of the operation of two turbines of the same characteristics; 
40 FIG. 23 is a characteristic showing a relation of a major steam flow amount at the operation of the turbines to a gen- 
eration output; 

FIG. 24 is a view showing a concept of a local optimal solution and global optical solution to the non-convexity of a 
valve point characteristic having turbines; 

FIG. 25 is a view showing a relation of a linear solution, non-linear solution and true optimal solution to the non- 
45 convexity of a valve point characteristic in turbines; and 

FIG. 26 is a block diagram showing a schematic form of an optimal operation system of a conventional power gen- 
eration plant. 

[Best Mode of Carrying Out the Invention] 

50 

A first embodiment of the present invention will be explained below with reference to the accompanying drawings. 
First, a brief explanation will be given of the present invention. 

As a plant system to which a plant utility optimization system is applied, an explanation will be given below of a non- 
utility generation system. The non-utility generation system is adapted to reutilize steam discharged from turbines and, 
55 there, a demand is made for both a steam discharged from the turbines and an energy output from the generators. 

The plant utility optimizing system of the present embodiment is directed to, while satisfying the steam demand and 
power demand as well as various kinds of restriction requirements, satisfying an operation state at a minimal total cost 
through the utilization of a dynamic programming, that is, in such a way as to allow the minimization of the total cost 
required for the operation of the plant system. 
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The above-mentioned dynamic programming is one procedure for solving a decision problem in which an optimal 
solution is normally decided in a multistage fashion with a result of its analysis providing a precondition to its sequential 
decision problem. In the dynamic programming, a decision selection sequence for a realizable optimal solution is called 
a "policy" and the policy of maximizing a specific function (criterion function) of a final state variable is called as an "opti- 
mal policy". The dynamic programming utilizes the "principle of optimality" throughout in which the optimal policy has 
the property that, whatever an initial state and initial decision are, a subsequent decision sequence has to be given the 
optimal policy with respect to a state produced as a result of the initial decision. The dynamic programming is set out in 
detail, for example, in TDynamic Programming, Princeton University Press, 1957J by R. Bellman. 

In a real plant dealt with in the present embodiment it is necessary that an optimal solution be found under the com- 
plex restriction conditions because of the ^utilization of the steam, the presence of a boiler group, the presence of var- 
ious steam systems, etc. It has not been possible to derive an optimal solution by the dynamic programming thus far 
employed. 

With the plant utility optimizing system of the present embodiment, not only means for optimizing the output of the 
turbine group but also various optimizing means are sequentially combined in an organic way along the arrangement 
of the plant and, by doing so, it is possible to find an optimal operation state of the plant satisfying a conplex demand. 
Briefly stated, first the processing utilizing the dynamic programming is performed twice and processing utilizing the lin- 
ear programming is performed once to find a solution (that is, an optimal solution candidate) for optimizing the output 
of the turbine group corresponding to respective operation state and then decide, as an optimal solution, an operation 
state for minimizing the total cost. 

Now the plant utility optimizing system of the present embodiment will be explained below in more detail. 

FIG. 1 is a diagrammatic view showing a plant energy system according to the present embodiment As shown in 
FIG. 1 , the system comprises a plant system including non-utility generation system 1 1 of a BTG (a boiler/turbine/ gen- 
erator), an external power input apparatus 12 and steam system 13, a plant controller 4 for controlling the non-utility 
generation apparatus 1 1 and external power input apparatus 12, and a plant utility optimizing system 2 connected to 
the plant controller 4 to find an optimal operation point (distribution amount) of the respective components in the plant 
system 6 and transfer the optimal operation distribution amount to the plant controller 4. 

In place of being connected to the plant system 6. the plant utility optimizing system 2 may be of such a type that 
the inputting and outputting of the data to and from the plant utility optimizing system 2 are done through the use of an 
information memory medium such as a magnetic disk. 

FIG. 2 diagrammatically shows one example of the flow of processing regarding the plant utility optimizing system 

2. 

Detailed system configuration information of the plant system 6 is input to the plant utility optimizing system 2 (step 

Step S12 decides a power demand/steam demand and inputs its value to the plant controller 4. 

Step S13 transmits the power demand/steam demand from the plant controller 4 to the plant utility optimizing sys- 
tem 2 and makes a request to execute optimization processing. 

The plant utility optimizing system 2 searches an optimal solution on the given power/steam demand (step S14). 

The plant utility optimizing system 2 transmits the searched optimal solution to the plant controller 4 (step S15). 

Based on the optimal solution given from the plant utility optimizing system 2 the plant controller 4 controls the plant 
system (step Si 6). 

FIG. 3 shows a detail of the plant system 6 in FIG. 1 . 

The present embodiment is equipped with four steam generation boilers (hereinafter referred to as boilers) B1 to 
B4, five steam turbines (hereinafter referred to as turbines) T1 to T5 and five generators G. 

The respective boilers B1 to B4 generate steam with the use of fuel fed. Let it be supposed that fuel amounts fed 
to the respective boilers B1 to B4 are given by f1 to f4 and that steam amounts generated by the boilers B1 to B4 are 
given by S1 to S4. 

L1 and L2 show high-pressure steam headers of different pressures (L1 is on the high voltage side). Further, S5 
and S6 show an input steam amount and output steam amount of a pressure reducing valve 20. 

The respective turbines T1 to T5 have functions p1 , p2, p3 and p5, respectively, of an input steam amount-output 
power characteristic and generate electric powers corresponding to the input steam amounts. The steam amounts input 
to the respective turbines T1 to T2 are represented by S1 and S5. 

The steams discharged from the respective turbines T1 to T5 are connected to corresponding load steam headers 
L3, L4 of different pressures. And a steam amount S1 ' is output from the steam system L3 and a steam amount S2 ? from 
the steam system L4. 

The electric powers generated from the generators G driven by the corresponding turbines T1 to T5 are output as 
a total power. 

The fuel amounts f 1 to f4 input to the corresponding boilers B1 to B4 have their upper and lower limit values deter- 
mined, respectively. Incidentally, further restriction is sometimes added as upper and lower limit values of those 
amounts of fuel input based on the amounts of fuels currently burned at the respective boilers B1 to B4. 
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Regarding the steam amounts S1 to S5 input to the turbines T1 to T5 their upper and lower limit values are deter- 
mined. 

In such a plant system 6, steam demands SD1 , SD2 (a minimal value to be satisfied by S1\ a minimal value to be 
satisfied by S2') for the respective steam systems L3, L4, as well as an energy demand ED, are determined. In the case 
s where the energy demand is not satisfied with the total generation power, an external commercial electric power is input 
for replenishment. 

In the case where the total generation amount of all the generators G exceeds the energy demand, a measure is 
taken to discard any excess portion or sell the generated power. 

Although in the plant system 6 of the present embodiment the turbines T1 to T5 have been explained as being not 

10 of a bleeding type, a plant system regarding the bleeding turbines can be so configured as to incorporate the bleeding 
turbines into the same arrangement as shown in FIG. 3 where these bleeding turbines are regarded as being connected 
in parallel with a plurality of turbines of non-bleeding type. 

FIG. 4 shows a schematic arrangement of a major portion of the plant utility optimizing system 2 according to the 
first embodiment of the present embodiment. As shown in FIG. 4, the plant utility optimizing system of the present 

15 embodiment comprises a fluid demand system optimizing section 21 , list converting section 22, list holding section 23, 
fluid input system optimizing section 24, energy/fluid demand input section 25, plant data holding section 26, fluid gen- 
eration optimizing section 27, list retrieving section 28 and a optimal solution outputting section 29. 

The plant data holding section 26 stores, therein, the upper and lower limit values of the fuel amounts f 1 to f4 input 
to the respective boilers B1 to B4, upper and lower limit values of steam amounts (steam discharge amounts) input to 

20 the respective turbines T1 to T5, functions p1 , p2, p3, p4 and p5 of an input steam amount-output power characteristic, 
and other various constants and functions such as functions of the input/output characteristic of the respective boilers. 
The plant data holding section 26 supplies, to the fluid demand system optimizing section 21 , list conversion section 22, 
fluid input system optimizing section 24 and fluid generation system optimizing section 27, their demanded numerical 
values and functions as will be set out below. 

25 First, the fluid demand system optimizing section 21 decides steam discharge amounts S1 , S2 and S4 of the 
respective turbines T1 . T2 and T4 through the use of the dynamic programming to maximize the total output power of 
the generators G driven by the turbines T1 , T2 and T4 connected to the steam system L3 and, in this case, the steam 
flow amounts in the steam system L3 are given their respective values. 

FIGS. 5 and 6 show one example of the operation of such fluid demand system maximizing section 21 . In FIGS. 5 

30 and 6, p1 , p2, .... p5 show the characteristic functions of the respective turbines and M1 , M2 and M4 show the upper 
limit values of the discharge steam amounts of the respective turbines T1 , T2, T4. For brevity in explanation it is 
assumed that M 1 , M2 and M4 are integers and the lower limit values of the steam amounts discharged from the respec- 
tive turbines T1 , T2, T4 all zeroes. Further, T and T represent arrays, a first [ ] following the T and T represents a row 
and the next [ ] represents a column, the following [0] [1] [2] representing first, second and third elements in the row/col- 

35 umn. 

Here, the steam flow amounts are found for the case where they are given respective integer values from 0 to M1 
+ M2 + M4. 

A table as shown in FIG. 7 is obtained by the operation as shown in FIGS. 5 and 6. 

Similarly, for the case where the steam flow amounts of steam systems L4 are given respective values, the fluid 
40 demand system optimizing section 21 decides steam discharge amounts S3, S5 of the respective turbines T3, T5 
through the use of the dynamic programming to maximize the total output power of the generators G driven by the tur- 
bines T3, T5 connected to the steam system L4. 

The operation of such fluid demand system optimizing section 21 is such as to correct M1 and M2 in FIG. 5 to M3 
and M5 and correct I" to Aj to "an end" in the case of NO in step S27. 
45 Even in this, the steam flow amount is found in the case where it is given a respective value from 0 to M3 + M5. 

For the respective steam flow amounts of the steam system L4, through the operation of the fluid demand system 
optimizing section 21, a table as shown in FIG. 8 is obtained where the optimal steam distribution and output power of 
the turbines T3, T5 are obtained. 

The tables shown in FIGS. 7 and 8 are sent to the next list conversion section 22. 
so Here, given that the steam amounts supplied from the boilers B1 to B4 to the L1 and L2 systems are <r1 , o2, then 

a1 = S1 + S2 

a2 = S3 + S4. 

55 

Thus, the following relations are established. 

S1 + S2 + S3 = <j1 - S5 
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S4+S5 = a2 + S6 
Here a linear relation, S5 = X • S6, exists between S5 and S6, an equation 

5 U + V = O 

is established. 
Here, 

10 u = S1 + S2 + X • S4 

v = S3 + A-S5 
a = a1 + X • a2 

where o represents the total amount of the generated steam in terms of the high pressure steam. 
The list conversion section 22 prepares a table as shown in FIG. 9 with the value u added to each column of the 
table shown in FIG. 7 and a table as shown in FIG. 10 with the valve v to the table as shown in FIG. 8. The value uM in 
FIG. 9 shows a maximal step value not exceeding M1 + M2 + X • M4 , while the value vM in FIG. 10 shows a maximal 
20 step value not exceeding M3 + X • MS . The list conversion section 22 sends these Tables as shown in FIGS. 9 and 10 
to the list holding section 23. 

Incidentally, simply preparing the tables as shown in FIGS. 9 and 10 by adding the values u and v calculated from 
the respective steam amounts in the steam system L3 and respective steam amounts in the steam system L4 to the 
tables as shown in FIGS. 7 and 8 results not necessarily in the values u and v being arranged at regular intervals. 
25 It is preferable that, in order to make the values u and v a constant step width, the list conversion section 22 decide 
the step widths of the values u and v and properly interpolate the steam amount, optimal steam distribution and output 
voltage with respect to the value u and v on the basis of the values of FIGS. 7 and 8. As the interpolation method the 
following may be conceived for instance. If the total steam amount corresponding to the values u and v exists on the 
row of the steam amount of L3, L4 on the tables shown in FIGS. 7 and 8, a corresponding column is adopted as it is 
30 and, if such is not the case, a row of a minimal value of those equal to, or exceeding, the total steam amount is adopted. 
This is another linear interpolation method by which, if the total steam value corresponding to the values and v 
exists on the table of FIGS. 7 and 8, a corresponding column is adopted and, if such is not the case, use is made of a 
column of a minimal value of those exceeding the total steam amount and column of a maximal value of those less than 
the total steam amount. 

35 The list holding section 23 stores, therein, the tables shown in FIGS. 9 and 10 on a given form. 
At the time when the plant system is decided, processing up to this can be performed. 

The energy demand and steam demand made from the plant controller 4 is sent from the energy/liquid demand 
input section 25 to the liquid input system optimizing section 24. Upon receipt of the already prepared tables as shown 
in FIGS. 9 and 10, as set out above, from the list holding section 23, the fluid input system optimizing section 24 
40 extracts, from the table as shown in FIG. 9, the row in which the steam amount of the steam system L3 satisfies the 
steam demand SD1 tor the L3 system and, similarly, from the table as shown in FIG. 10, the column in which the steam 
demand SD2 satisfies the steam demand SD2 for the L4 system. 

By this operation it is possible to obtain the table as shown in FIGS. 1 1 and 12. It is to be noted that, in the table as 
shown in FIG. 1 1 , L3I and L3u are lower and upper limit values of the L3 steam amount as shown in FIG. 9 they satisfy 
45 the steam demand for the L3 system and Urn, (sm1 , sm2, sm4) and Em3 correspond to u, (sM1 , sM2, sM4) and output 
power in the column of L3I. The same thing is also true of uM (sM2, sM2, sM3). The meaning of the respective symbols 
in Table as shown in FIG. 12 is the same as shown in Table in FIG. 1 1 . 

For the case of letting the total input steam amount be respective values in terms of a high pressure level, the fluid 
input system optimizing section 24 uses the dynamic programming based on the tables as shown in FIGS. 11 and 12 
so to decide steam supply amounts S1 to S5 to the respective turbines T1 to T5 so that the total output power of the gen- 
erators G driven by the turbines T1 to T5 may be made maximal. 

The tables of FIGS. 1 1 and 1 2 are regarded as two-dimensional arrays, 4x11,4x12, and referred to as A, B, not- 
ing that 11,12 represent the lengths of the tables. For example, A [0] [1 1 -1 2] corresponds to uM. 

FIG. 13 shows the operation of the fluid input system optimizing section 24. In FIG. 13. T represents an array in 
55 which a first [ ] and next [ ] represent the row and column, respectively, with the following [0], [1] representing first and 
second elements, respectively, in the row/column array. 

An A [2] D] + B [2] [k] in FIG. 13 shows a merge operation of a list A [2] Q] and let B [2] [kj. For example, let the 
elements in the list A [2] [j] and these in the list B [2] [k] be (a1 , a2, a4) and (a3, a5), respectively, and as a result of the 
merge operation these elements of A [2] [j] + B [2] [k] becomes (a1 , a2, a3, a4, a5). 



10 

BNSDOCtD: <EP 0731397A1_I_> 



EP 0 731 397 A1 



By the operation as shown in FIG. 13 It is possible to obtain, as an array T r a Table as shown in FIG. 14. This table 
is sent to the fluid generation optimizing section 27. 

Here, in the table as shown in FIG. 14, the steam distribution (s1, s5) with the total power output maximal to 
respective column (that is, the respective input steam amount) is obtained, so that the steam amounts a1 , o2 flowing 
s through the respective systems of L1 , L2 with respect to the respective total input steam amounts are uniquely decided. 

Under the condition that the lower and upper values of amounts of fuel, f1 to f4, are satisfied with respect to the 
respective columns in the table as shown in FIG. 14, that is, with respect to the respective total input steam amounts 
and that the steams of a1 , o2 are supplied to the respective systems of L1 , L2, the fluid generation system optimizing 
section 27 finds a fuel distribution by which the fuel cost is minimized. The boiler characteristics of the respective boilers 
io B1 to B4 become, in any case, a linear function of the fuel costs also become a linear function of the fuel. Therefore, 
the fuel distribution and fuel cost can be found by the linear or the nonlinear programming. 

By doing so, those values SI to S6 and f 1 to f4 are found with respect to the respective columns of the table as 
shown in FIG. 14. The fluid generation system optimizing section 27 prepares a table as shown in FIG. 15 while adding 
these values to the table as shown in FIG. 14. 
15 The table is sent to the list retrieving section 28. 

At this time, it is possible to obtain operation parameters with which the total power output is maximal with respect 
to the respective fuel cost. That is, the respective columns can be regarded as candidates for an optimal solution to be 
sought. 

The list retrieving section 28 decides an optimal solution from the respective optimal candidates of the table as 
20 shown in FIG. 15. 

First, the list retrieving section 28 calculates the total costs, for each column, that is, for each optimal solution can- 
didate, on the basis of the values on the Table as shown in FIG. 15. 
The total cost is determined, for example, by the following equation. 

25 (total cost) = (total fuel cost) + e (ED - (total output power)) 

where 

e: unit price of the electric power to be bought 
30 ED: energy demand. 

This case is based on the condition that the total output power amount is set to be lower than the energy amount. 
In the case where the total output power amount exceeds an energy demand, it may set, for example, (ED - (total output 
power)) to be zero. 

35 Or in the case where any excess portion over the energy demand can be sold, e' may be calculated as a unit value 
to be sold. 

In the case where ED a= the total output power, it may be calculated as 

(the total cost) = (total fuel cost) + e (ED - (total output power)). 

40 

In the case where ED < total output power, it may be calculated as 

(total cost) = (total fuel cost) + e* (ED - (total output power)). 

45 The list retrieving section 28 prepares a table as shown in FIG. 16 while adding the calculated total cost to the 
respective column of the table as shown in FIG. 15. From this table a column of a minimal total cost is found and a set 
of various kinds of values in the column is supplied as an optimal solution on the whole system to the optimal output 
section 29. 

The optimal solution output section 29 holds the optimal solution sent from the list retrieving section 28. Out of var- 
50 ious values in the optimal solution, requisite information is sent as an optimal operation parameter to the plant controller 

4. 

According to the plant utility optimizing system of the first embodiment, it is possible to satisfy various demands for 
the plant system and find an optimal operation state reflecting the non-convex characteristic of the turbines. Therefore, 
the operation cost can be largely reduced in comparison with the conventional counterpart. 
55 Although the first embodiment has been explained in connection with the operation of the plant utility optimizing 
system fit for the plant system of FIG. 3, the present invention can also be applied to the case where the number of the 
boilers and turbines as well as the kinds of steam systems are changed. 

In the case where the electric power and steam can be supplied to the power load and steam load by the non-utility 
generation equipment without the need to buy electric power from an outer power source, it may be possible to, instead 
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of performing the processing by the list retrieving section 28, find an optimal solution (optimal operation state taking the 
non-convex characteristic into consideration) on the basis of the table prepared by the fluid generation system optimiz- 
ing section 32. 

Further, even if other elements are added to the plant system as shown in FIG. 3, an optimal solution can be 
obtained through the proper correction of the plant utility optimizing system. The present invention can be applied not 
only to the non-utility generation system but also generally to the plant system including those elements having the non- 
convex characteristic. 

A second embodiment of the present invention will be explained below. 

FIG. 17 is a block diagram showing a schematic form of an optimal operation system of a power generation plant 
(BTG plant). 

In FIG. 1 7, 26 shows a plant model separated into the turbine-and the boiler-side model. In the plant model 26, the 
turbine characteristic, turbine operation restrictions, steam balance, etc., are given to the turbine-side model and the 
boiler operation, boiler operation restrictions, steam balance, etc., are given to the boiler-side model. 

21 shows a turbine output maximizing point calculation section for deciding the steam distribution of each turbine 
to maximize the total of turbine outputs on the basis of data of each steam amount supplied from the turbine-side model 
and 22 shows a fuel cost minimizing point calculation section for deciding a steam distribution amount of each boiler to 
generate the total of per-pressure supply steam amounts at a minimal cost the supply steam amounts being given from 
the boiler-side model at each point found at the turbine output maximizing point calculation point 21 . 

Further, 23 shows a total energy cost minimizing point calculation section for searching a total energy cost minimiz- 
ing point based on the steam distribution amount of each boiler, that is, on the steam distrtoution amount given as a 
plant load (electric load amount, steam load amount) of a plant state, including a power buying unit price, and deter- 
mined by the fuel cost minimizing point calculation section 23. 

On the other hand, 24 shows an editor for changing the plant model 26. The editor 24 is comprised of a plant model 
preparing means for preparing/editing system graphics for modifying or correcting the plant model used at the calcula- 
tions of the turbine output maximizing point calculation section 21 and fuel cost minimizing point calculation section 21 
as well as the system model described with mathematical expressions. 25 shows a data generator for enabling data 
which is necessary for the turbine output maximizing point calculation section 21 and fuel cost minimizing point calcu- 
lation section 22 to be changed from an optimalized model modified or corrected by the editor 24 to a format matched 
to the calculation system. 

The function of the maximal operation system of the BTG plant thus arranged will be explained below. 

The turbine output maximizing point calculation section 21 decides the steam distribution of each turbine to maxi- 
mize a turbine output at each supply steam amount total given from the turbine side model. In this case, in the plant 
including a combination set of turbines having a non-convex valve point characteristic, the supply steam amount total 
is imparted to the turbines and a steam distribution combination is found at that turbine whose power generation 
amount is maximal at that time. 

Here, in the case of finding the steam distribution combination whose power generation amount is maximal relative 
to the turbine, the most primitive method for finding a global optimal point is by dividing the supply steam amount total 
into a few hundred steps, varying the steam amounts distributed to the respective turbines at each one step to find 
power generation output totals and selecting, out of these, a distribution combination whose power generation is max- 
imal, noting that, in order to prevent a lowering in calculation accuracy, the number of steps in the step function for one 
turbine is adequately greater (for example, 1 00) than about 3 to 5 for a turbine case (the number of valve points < indus- 
trial steam turbines). In such all point checking, a vast number of times of calculations is necessary and a difficulty is 
encountered in actual practice. 

In this situation, a DP (dynamic programming: subsequent optimal method) is used, according to the present inven- 
tion, as a means, for example, for obtaining a maximal point in less number of times, a method which is the same as an 
"all points" trial. 

The result of calculation by the DP is found, for each a few hundreds to a few thousands of input steam amounts 
corresponding to sums of step functions of all the turbines over a range from a minimal value to a maximal value taken 
by the sums of all the turbines for example. 

Further, the fuel cost minimizing point calculation section 22 decides a steam distribution amount of each boiler so 
as to generate, at a minimal cost, the sum total of per-pressure supply steam amounts given from the boiler-side model 
at each point of a turbine input steam amount sum total found at the turbine output maximizing point calculation section 
21 . In this case, the fuel unit is not varied in a time unit as in the case of an electric power unit and is given a fixed value 
and. if the need arises, can be varied with a manual instruction. 

Since the boiler characteristic is substantially flat as set out above, calculation is made at the fuel cost minimizing 
point calculation section 22, through the use of a linear approximation, with less affected accuracy in a resultant solu- 
tion, so that it is Possible to apply an LP which is quicker in calculation speed. 

With the use of a fuel cost minimum/power generation amount maximum distribution table at each point of the tur- 
bine input steam amount total found at the fuel cost minimizing point calculation section 22, the total energy minimizing 
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point calculation section 23 finds a total energy cost, at the time of giving a load power amount, with a per steam amount 
sum total as follows: 

C = Et (PI - Pg) + Cf 

5 





wherein 


C: 


total energy cost (yen/H) 


Et: 


buying power unit (yen/KWH) 


PI: 


load power (KW) 


Pg: 


generation power (KW) 


Cf: 


fuel cost (yen/H) 



An optimal distribution point is found by rearranging the costs C of the total energies in a decreasing value order 
15 and, checking to see whether or not the load supply steam amount sum total of each pressure satisfies the load steam 
amount sum total and finding such an initially satisfied point as that optimal distribution point. 

By the operation as set out above, it is possible to obtain a power/steam optimal distribution point with respect to 
the power load and steam load per pressure. 

In the case where the optimizing model is corrected with respect to the system change and apparatus's expansion 
20 of the target energy plant as well as the apparatus specification, characteristic variation, etc., the user prepares a steam 
system diagram in the editor 24 in an "AD" feeling and inputs the apparatus characteristic, rating, etc., from the steam 
system diagram to correct the optimizing model. 

In this case, it is possible for the user to edit the operation restriction condition, equality restriction, inequality- 
restricted conditional equation through the editor in the "word processing" feeling. 
25 Upon receipt of the system diagram, apparatus data and restricted conditional equation prepared by the editor 24, 
the data generator 25 automatically generates data necessary for optimization calculation from these information items 
in a form fit for an solution, so that the plant model 26 is corrected. 

In this case, the optimal calculation logic is separated from a data description portion such as the system configu- 
ration, characteristics and restriction condition and, with the use of a system not requiring any timing parameters in the 
30 optimal solution method, the system maintenance can be achieved simply by modifying the model. 

By separating the steam distribution calculating means relative to the turbine having the non-convex characteristic 
and generation steam distribution calculation means of the boiler having the convex characteristic and allowing ready 
linear approximation it is possible, according to the present invention, to, upon calculation of the steam distribution to 
the turbine, apply the sequential calculation procedure (DP), allowing a global optimal solution equivalent to an all-point 
35 search, which is carried out with far less number of times than that necessary for the all point search. It is, therefore, 
possible to achieve an improved accuracy solution and, further, to employ a high-speed linear programming (LP) rela- 
tive to the generation steam distribution of the boiler. 

By adopting such a two-step processing it is possible to correctly get the global optimal point including a non-con- 
vex turbine valve characteristic and, in addition, achieve a system not requiring any special search logic and parameter 
40 adjustment, so that optimal calculation can be carried out with the use of the user's corrected model. 

According to the second embodiment of the present invention a solution can be obtained with higher accuracy than 
according to the conventional system. Further, it is not necessary to tune the calculation parameters. Thus the user of 
no expert's knowledge on the optimization calculation can carry out the system modification of the target energy plant 
as well as the system maintenance such as the correction of the plant model corresponding to the apparatus modif ica- 
45 tion. 

Claims 

1. A utility optimizing method for searching optimal operation parameters of a plant including a plurality of first plant 
so elements having a linear input/output characteristic with first physical amounts as an input and second physical 
amounts as an output and a plurality of second plant elements having an input/output characteristic with the second 
physical amounts as an input, corresponding to the output of the first plant element, and third physical amounts as 
an output, the input/output characteristic containing a non-convex characteristic portion, the method comprising: 

55 a first processing step of finding a distribution of the second physical amounts to the second plant elements by 

a dynamic programming so as to maximize the sum total of the third physical amounts corresponding to the 
sum total of the second physical amounts in a predetermined range satisfying a demand of the second physical 
amounts and output from the second plant elements; 
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a second processing step of finding a distribution of the second physical amounts to the second plant elements 
by the dynamic programming on the basis of the sum total of the second physical amounts, sum total of the 
third physical amounts and distribution of the second physical amounts to the second plant elements obtained 
by the first processing step to maximize the sum total of the third physical amounts corresponding to the sum 
5 total of the second physical amounts to the second plant elements and output from the respective second plant 

elements; and 

a third processing step of finding, by a linear programming or nonlinear programming, the first physical 
amounts input to the second plant elements and evaluation values, necessary for supplying the first physical 
amounts, with respect to the sum total of the second physical amounts, sum total of the third physical amounts 
10 and distribution of the second physical amounts to the second plant elements obtained at the second process- 

ing step. 

2. A utility optimizing method for searching optimal operation parameters of a plant including a plurality of first plant 
elements having a linear input/output characteristic with first physical amounts as an input and second physical 
15 amounts as an output and a plurality of second plant elements having an input/output characteristic with the second 
physical amounts as an input, corresponding to the output of the first plant elements, and third physical amounts 
as an output, the input/output characteristic containing a non-convex characteristic portion, the method comprising: 

a first processing step of finding a distribution of the second physical amounts to the second plant elements by 
20 a dynamic programming to maximize the sum total of the third physical amounts corresponding to the sum total 

of the second physical amounts in a predetermined range satisfying a first demand for a group of such second 
physical amounts and output from the second plant elements associated with the group; 
a second processing step of finding a distribution of the second physical amounts to the second plant elements 
by a dynamic programming on the basis of a set of the sum total of the second physical amounts, sum total of 
25 the t hir d physical amounts and distribution of the second physical amounts to the second plant elements 

obtained by the first processing step at each second physical amount in a group to maximize the sum total of 
the third physical amounts corresponding to the sum total of the second physical amounts to all the second 
plant elements and output from said all the second plant elements; 

a third processing step of finding the first physical amounts input to the respective second plant elements and 
30 a first cost required to supply the first physical amounts, by a linear programming or nonlinear programming, 

with respect to each set of the sum total of the second physical amounts, sum total of the third physical 
amounts and distribution of the second physical amounts to the second plant elements obtained by the second 
processing step; and 

a fourth processing step of, on the basis of each set of the sum total of the third physical amounts and first cost 
35 obtained at the third processing step, calculating a second cost required to receive a power which corresponds 

to a shortage in a second predetermined power demand for the total sum of the third physical amounts from 
an outside source to find a total cost for each set and deciding, as an optimal solution, each value belonging 
to a set in which the total cost is minimized. 



40 3. The method according to claim 2, wherein when, in the fourth processing step, the sum total of the third physical 
amounts exceeds the second power demand, the second cost is made zero or the second cost is calculated as a 
negative value when the second cost is obtained by supplying an excess of the sum total of the third physical 
amounts over the second power demand to the outside. 



45 4. A utility optimizing method for searching optimal operation parameters of a power generation plant including a plu- 
rality of boilers, a plurality of turbines operated by steams generated from the boilers, and generators driven by the 
respective turbines to generate electric power in which the electric power is supplied to electric power loads and 
the steams bled from the respective turbines are supplied to steam loads via a plurality of turbine discharge sys- 
tems, the method comprising: 

50 

a first processing step of finding a steam flow amount distribution to the respective turbine by a dynamic pro- 
gramming to maximize the sum total of generation output powers of the turbines corresponding to respective 
steam flow amounts in a predetermined range satisfying a steam demand to the respective turbine discharge 
system and connected to the same turbine discharge system; 
55 a second processing step of, on the basis of a set of the steam flow amount, generation output power and 

steam flow amount distribution obtained by the first processing step at each turbine discharge system, finding 
a steam flow amount distribution to the respective turbine by a dynamic programming to maximize a total gen- 
eration output power of all the turbines corresponding to the total input steam flow amounts to all the turbines; 
and 
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a third processing step of finding a supply fuel distribution to the respective boilers and total fuel cost by a linear 
or a nonlinear programming with respect to each set of the total input steam flow amount, total generation out- 
put power and steam flow amount distribution to the respective turbines. 

5 5. A utility optimizing method by searching optimal operation parameters of a power generation plant including a plu- 
rality of boilers, a plurality of turbines operated by steams generated from the boilers and generators driven by the 
respective turbines to generate electric power in which the electric power is supplied to electric power loads while 
the steams bled from the respective turbines are supplied to steam loads via a plurality of turbine discharge sys- 
tems, 

10 

a first processing step of finding a steam flow amount distribution to the respective turbine by a dynamic pro- 
gramming to maximize the sum total of the generation output powers to the turbines corresponding to the 
steam flow amounts in a predetermined range satisfying a steam demand to the respective turbine discharge 
system and connected to the same turbine discharge system; 

75 a second processing step of finding a steam flow amount distribution to the respective turbines by the dynamic 

programming on the basis of a set of the steam flow amount, generation output power and steam flow amount 
distribution obtained by the first processing step at each turbine discharge system to maximize a total genera- 
tion output power of the turbines corresponding to a total input steam flow amount to all the turbines; 
a third processing step of finding, by a linear or a nonlinear programming, a supply fuel distribution to the 

20 respective boilers and total fuel cost with respect to each set of the total input steam flow amount, total gener- 

ation output power and steam flow amount distribution to the respective turbines; and 

a fourth processing step of, on the basis of each set of a total generation output power and total fuel cost 
obtained by the third processing step, calculating an external cost required to receive a power which corre- 
sponds to a shortage in a predetermined power demand for the total generation output power from an outside 
25 source to find a total cost for each set and deciding, as an optimal solution, each value belonging to a set in 

which the total cost is minimized. 

6. The method according to claim 5 wherein when, in the fourth processing step, the total generation output power 
exceeds the output demand the external cost is made zero or the external cost is calculated as a negative value 

30 when the second cost is obtained by supplying an excess of the total generation output power over the power 
demand an outside. 

7. In a power generation plant including a plurality of boilers, a plurality of turbines operated by steams generated from 
the boilers and generators driven by the respective turbines to generate electric power wherein the electric power 

35 is supplied to electric power loads while the steams fled from the respective turbines are supplied to the steam 
loads, a plant utility optimizing system comprising: 

a plant model prepared by separation into a turbine-side and a boiler-side model; 

turbine output maximum point calculating means for deciding such a steam distribution of the turbines as to 
40 maximize a total of outputs of the respective turbines for each supply steam amount given from the turbine-side 

model in the plant model; 

fuel cost minimizing point calculating means for deciding such a steam distribution amount of each boiler as to 
generate at a minimal cost, the sum total of per pressure supply steam amounts given from the boiler-side 
model at each point found by the turbine output maximum point calculating means; and 
45 total energy cost minimizing point calculating means for searching a total energy cost minimum operation point 

to electric loads and steam loads on the basis of a steam distribution amount decided by the fuel cost minimiz- 
ing point calculating means. 

8. In a power generation plant including a plurality of boilers, a plurality of turbines operated by steams generated from 
so the boilers and generators driven by the respective turbines to generate electric power in which the electric power 

is supplied to electric power loads while steam bled from the respective turbine is supplied to steam loads, a plant 
utility optimizing system comprising: 

a plant model prepared by separation into a turbine-side and a boiler-side model; 
55 turbine output maximum point calculating means for deciding such a steam distribution of the turbines as to 

maximize a total of outputs of the respective turbines for each supply amount given from the turbine-side 
model; 
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fuel cost minimizing point calculating means for deciding such a steam distribution amount of the respective 
boilers as to generate, at a minimal cost the sum total of per pressure supply steam amounts given from the 
boiler-side model at each point found by the turbine output maximum point calculating means; 
total energy cost minimizing point calculating means for searching a total energy cost minimum operation point 
to electric power loads and steam loads on the basis of a steam distribution amount of the respective boilers 
decided by the fuel cost minimizing point calculating means; 

plant model preparing means for preparing/editing system graphics for modifying or correcting plant models 
used under the calculation of the tuibine output maximum point calculating means and fuel cost minimizing 
point calculating means and system models described with a mathematical expression; and 
a data generator for generating data necessary for optimal calculation from the model prepared by the plant 
model preparing means and for giving the data to the turtine-and boiler-side models in the plant model. 
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PLANT CONTROLLER 



NON- UTILITY GENERATION 
SYSTEM (BTG) 



EXTERNAL 
POWER INPUT 



1 2 



1 1 



STEAM SYSTEM 



1 3 



PLANT UTILITY OPTIMIZING SYSTEM 



FIG. \ 



6' 



1 



S 1 1 



INPUT PLANT CONFIGNRATION TO 

A PLANT UTILITY OPTIMIZING SYSTEM 



1 * S12 

DESIDE POWER AND STEAM DEMANDS INPUT THE POWER 
AND STEAM DEMANDS INTO THE CONTROLLER 



1 3 



> 

TRANSMIT THE POWER ANC 
THE CONTROLLER TO AN OF 
MAKE AN OPTIMIZATION RE 


) STEAM DEMANDS FROM 
>TIMIZING SYSTEM-, 
QUEST FROM CONTROLLER 


\ 





SEARCH AN OPTIMAL SOLUTION 



1 5 



> 

TRANSMIT THE OPTIM 
OPTIMIZING SYSTEM T 


AL SOLUTION FROM AN 
0 A CONTROLLER 




<S16 


CONTROLLER CONTROLS AN ENERGY 
SYSTEM OF THE PLANT 



F I 6. 2 



17 



BNSDOCID: <EP 0731397A1J_> 



EP 0 731 397 A1 




BNSOOCID: <EP 0731397A1 J_> 



EP 0 731 397 A1 



ELECTRIC POWER/ 
FLUID DEMAND 
INPUT SECTION 



2 5 



T 



2 1 



1 



FLNID DEMAND 
SYSTEM OPTIMIZING 
SECTION 



2 2 



1 



LIST CONVERSION 
SECTION 



23 



1 



LIST HOLDING 
SECTION 



FLUID INPUT SYSTEM 
OPTIMIZING SECTION 



2 4 



FLUID GENERATION 
SYSTEM OPTIMIZING 
SECTION 



27 



T 



LIST RETRIEVING 
SECTION 



2 8 



Y 



OPTIMAL SOLUTION 
OUTPUT SECTION 



2 9 



r 



PLANT DATA 
HOLDING SECTION 



T 



2 6 



F I 6. 4 



0731397A1 I 



19 



EP0 731 397 A1 



( START ) 




S23 





> 


YES 

f 


TCO)[j]= 


-10 0 0 0 \ 


S24x ; 


f 






j = j + 1 





1 



S2 5 



J =0 




1 



S2 9 



TC0][j+k] = p1(j)+p2(k) 
T[1][j+k][0]=j 
T[1][j+k][1] = k 



S32 



k = k+1 



FIG. 5 



20 



BNSDOCID: <EP 0731397A1_I_> 



EP0 731 397 A1 



S3 3- 




j=0 



S34< 



j £M1 +M2+M4? 



1 



S37 



S 3 5s 


YES 


j = 


= 0 




r co][j]=-i oooo 


> 






S 3 6s 




. <S38 


j = j + 1 


k = 0 






FIG. 6 



21 



EP 0 731 397 A1 



STEAM AMOUNT 
OF L3 


0 




Ml +M2+M4 


OPTIMAL STEAM 
DISTRIBUTION 
(SI.S2. S3) 


(0, 0, 0) 




(Ml. M2, M4) 


OUTPUT POWER 


0 




pi (Ml) + p 2 (M2) +p4 (M4) 



FIG. 7 



STEAM AMOUNT OF L4 


0 




M3+M5 


OPTIMAL STEAM 
DISTRIBUTION 
( S3 , s5 ) 


(0. 0) 


• • • 


(M3. M5) 


OUTPUT POWER 


0 




p3 (M3) + p 5 (M5) 



FIG. 8 



u 


0 




uM 


STEAM AMOUNT 
OF L3 


0 




M1 +M2+M4 


OPTIMAL STEAM 

DISTRIBUTION 

(S1.S2.S3) 


(0. 0. 0) 


• • • 


(Ml. M2. M4) 


OUTPUT POWER 


0 




p 1 (M1) +p2 (M2) +p4 (M4) 



FIG. 9 



V 


0 




vM 


STEAM AMOUNT OF L4 


0 




M3+M5 


OPTIMAL STEAM 
DISTRIBUTION 
( S3.S5) 


(0. 0) 




(M3, M5) 


OUTPUT POWER 


0 


• • • 


p3 (M3) +p5 (M5) 



FIG. \0 



BNSDOCID: <EP 0731397A1_L> 



22 



EP0 731 397 A1 



u 


Urn 




UM 


STEAM AMOUNT 
OF L3 


L3 I 




L3 u 


OPTIMAL STEAM 
DISTRIBUTION 
(SI.S2, S4) 


(smi, sm2, sm3) 




(sMl, sM2, sM4) 


OUTPUT POWER 


Em3 




EM3 



FIG. H 



V 


Vm 




VM 


STEAM AMOUNT OF L4 


L4 £ 




L4 u 


OPTIMAL STEAM DISTRIBUTION 
( Sl,s5) 


(sm3, sm5) 


• • • 


(sM3. sM5) 


OUTPUT POWER 


Em 4 




EM4 



F I G. 12 



TOTAL INPUT STEAM AMOUNT 
(IN TERMS OF HIGH PRESSURE ) 


Um+Vm 




UM+VM 


( STEAM OF L3, STEAM OF L4 ) 


(L3 L4 i) 




(L3u. L4u) 


OPTIMAL STEAM DISTRIBUTION 
(St, ,S5) 


(smi 

sm5) 




(sM 1 

sM5) 


TOTAL OUTPUT POWER 


Em3+Em4 


• • • 


EM3+EM4 



FIG. 14 



BNSDOCID: <EP 0731397A1_I_> 



23 



EP0 731 397 A1 



( START ) 




S5 5 



S5 3 



3 \ 




YES 


T[3HU= 


-10 0 0 0 


S54s , 


f 






j=j +1 





j = o 



1 



S5 6 



k = 0 



( END ) 




S6 0; 

AC3][j]+BC3][k]^TC3][j+k]? 



S 6 1 



1 



YES 



T[3][j+k]=AC3][j] + B[3]Ck] 
TC0][j+k]=AC0][j]+B[0]Ck] 
TCl3[j+k][0]=A[1][j] 
TCl][j+k][1] = BCl][k] 
T[2][j+k]=A[23[j]+B[2][k] 



S62 



1 



k = k+1 



FIG. 13 



24 



EP 0 731 397 A1 



TOTAL INPUT STEAM AMOUNT 
( IN TERMS OF HIGH PRESSURE ) 


Um+Vm 




UM+VM 


(STEAM OF L3, STEAM OF L4 ) 


(L3«, L4£) 




(L3iu L4u) 


STEAM DISTRIBUTION 
( SI, . . . , S5) 


(sml 

sm5) 




CsM1 

sM5) 


TOTAL OUTPUT POWER 


Em3+Em4 




EM3+EM4 


FUEL DISTRIBUTION 


(*, *, *, *) 




(*, *, *, *) 


TOTAL FUEL COST 


* 




* 



FIG. 15 



TOTAL INPUT STEAM AMOUNT 1 
(IN TERMS OF HIGH PRESSURE) 


Um+Vm 




UM+VM 


( STEAM OF L3 .STEAM OF L4 ) 


(L3 2. L4 1!) 




(L3u. L4u) 


STEAM DISTRIBUTION 
(SI,. . . ,S5 ) 


(sml 

sm5) 




(sMI 

sM5) 


TOTAL OUTPUT POWER 


Em3+Em4 




EM3+EM4 


FUEL DISTRIBUTION 


(*, *. *, *) 




(*. *. *, *) 


TOTAL FUEL COST 


* 




* 


TOTAL COST 


* 




* 



FIG. 16 



BNSDOCID: <EP 0731397A1_I_> 



25 



EP0 731 397 A1 




ro ^, 



o 

or o 
uj z 



LU M <Z 
h-COZ — JO 

po=0<uJ 

— * * 



10 r 



or 
o 

cr 



Z5-J0 



C\J 



cc 
o 

Q 
LJ 



o 
a.. 



L>=-IO 

UJZJU 
3-<UJ 



£ o 
o ^ 



o 
o cc 

— or — 
ft: i-O 

9?CK 2 
-J ^Ul< 
UJ 23 
Q OO < 



<D 



a. 



ODQQ< 

cr cr uj 



* 0 
1 cr I— 

1 LU CO 

!d E 

UJ 

_j 

K o 



2 o 
o «m 

O 

cr 
h- 

cn 

UJ 

cr 



o 

— uj 

£^ 
cr 3 

UJ 
O H- 
CD CO 




26 

BNSDOCID: <EP 0731397A1_I_> 



EP0 731 397 A1 




EP0 731 397 A1 




EP0 731 397 A1 




LOWER FUEL AMOUNT UPPER 
LIMIT (kl/H)— LIMIT 

FIG. 20 



UPPER 

STEAM AMOUNT (T/H ) —LIMIT 

FIG. Z\ 



FIG. 22 



FIG. 23 



50T/H 40T/H \ 



T 



50T/H =* 60T/H 



lOOT/H 



i=!L> 

55 



e> uj 



2I8 


2IO 








200 ~_ /j 










; 40 


1 50 


:60 



MAIN STEAM ( t/H)^ 
FLOW EFFICIENCY 



29 



BNSDOCID: <EP 0731397A1_I_> 



EP0 731 397 A1 




GLOBAL OPTIMAL POINT 
(TRUE OPTIMAL POINT) 



^ v y 

FIG. 24 



BNSDOCID: <EP 0731397A1 I > 



30 



EP0 731 397 A1 




F I 6. 25 



32 

cL 

PLANT MODEL 

TURBINE 
CHARACTERISTIC 

TURBIN OPERATION 
RESTRICTIONS 

BOILER 

CHARACTERISTIC 

BOILER OPERATION 
RESTRICTIONS 
STEAM BALANCE etc 



PLANT STATE 
ELECTRIC LOAD AMOUNT 
STEAM LOAD AMOUNT 
ENERGY BUYING UNIT PRICE 



31 

C±_ 

TOTAL ENERGY 
COST 

MINIMIZING POINT 

CALCULATION 
SECTION 

" TUNING- 1 
REQUIRED I 
CALCULATION j 
PARAMETERS 



OPTIMAL 
SOLUTION 



FIG. 26 



BNSDOCID: <EP 0731397A1J_> 



31 



EP0 731 397 A1 



INTERNATIONAL SEARCH REPORT 



International application No. 

PCT/JP95/01945 



A. CLASSIFICATION OF SUBJECT MATTER 
Int. CI* G05B13/02, F01K13/02 

According to International Patent Qasaifioabon (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 
Int. CI 6 G05B13/02, F01K13/02 



Documentation searched other than minimum documentation to the extent that such documents are ioduded in the fields searched 

Jitsuyo Shinan Koho 1926 - 19 95 

Kokai Jitsuyo Shinan Koho 1973. - 1 995 

Electronic dau base consulted during the international search (name of data base and, where practicable, search terms used) 



JOIS 



C DOCUMENTS CONSIDERED TO BE RELEVANT 



Category* 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



JP, 62-170709, A (Hitachi, Ltd.), 
July 27, 1987 (27. 07. 87) 

JP, 62-212701, A (Yokogawa Electric Corp. J, 
September 1.8, 19 87 (18. 09. 87) 

JP, 3-188562, A (Hitachi, Ltd.), 
August 16, 1991 (16. 08. 91) 

JP, 2-264303, A (Idemitsu Kosan Co., Ltd.), 
October 29, 1990 (29. 10. 90) 

JP, 2-49904, A (Mitsubishi Heavy Industries, 
Ltd . ) , 

February 20, 1990 (20. 02. 90) 

JP, 4-272407, A (Toshiba Corp.), 
September 29, 1992 (29. 02. 92) 

JP, 3-17204, U (Mitsubishi Heavy Industries, 
Ltd . ) , 



1 



8 
8 

1 - 8 
1 - 8 
I - 8 

1-8 
1-8 



[x] Further documents are listed in the continuation of Box C. Q See patent family annex. 



Spedal categories of died documents: 



later documeot published after the international Cling date or priority 
date and not to conflict with Use application but died to understand 
the principle or theory underlying the invention 



~ A " f °f»™«» defining the general state of the art which is not considered 

to be of particular rdevance . 

M E" earlier document but published on or after the international filing date **X M document of particular relevance; the cJ aimed ioveotioo cannot be 

**L" documeot which may throw doubts on orioritv cJaimfYl or which » considered novel or cannot be considered to involve an inventive 

Fredas P wb,ia,io0 °* <^<* 0 ofSS Wbea the doCUmenl tok *° « lo °* 

! Pe °* fe * S °r * pcoCed ) ~Y" document of particular relevance; the daimed invention cannot be 

u oocumcot referring to an oral disclosure, use, exhibition or other considered to involve an inventive step wbeo the document is 

lnc * ftS combined with one or more other such documents, such combination 

"P" document published prior to the international Wing dale but later than ***** obvious *° a P 61500 >ki,,ed ^ the art 

the pnonty dale daimed document member of the same patent family 



Date of the actual completion of the international search 

December 5, 1995 (05. )2. 95) 



Date of mailing of the international search report 

December. 26 , 1995 (26. 12, 



95) 



Name and mailing address of the ISA/ 

Japanese Patent Office 
Facsimile No. 



Authorized officer. 



Telephone No. 



Form PCT/ISA^10 (second sheet) (July 1992) 



32 

BNSOCCID: <EP 0731397A1_I_> 



# 



EP 0 731 397 A1 



# 



INTERNATIONAL SEARCH REPORT 



International application No. 

PCT/JP95/C1945 



C (Continuation). DOCUMENTS CONSIDERED TO BE RELEVANT 



Category* 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to daim No. 



July 3, 1989 (03. 07. 89) 

Youchi Kita "The Operation System for 
Optimizing BTG" Vol. 49, No. 1, (1995) 
Japan Association of the Pulp and Paper 
Industry, pp. 79-85 



1-8 



Form PCT/1SA/210 (continuation of second sheet) (July 1992) 



33 



BNSDOCID: <EP 0731397A1_I_> 



THIS PAGE BUm wm f 



♦ 



EP0 731 397 A1 



* 



m m j» m, m m m mm m ■ 



11 



PLANT CONTROLLER 



T 

Jl. 



HON -UTILITY GENERATION 
SYSTEM 1BTG1 



EXTERNAL 
POWER IMPUT 



\ 1 



1 1 



STEAM SYSTEM 



1 3 



PLANT UTILITY OPTIMIZING SYSTEM 



FIG. \ 



6' 



INPUT PLANT CONFIGURATION TO 

A PLANT UTILITY OPTIMIZING SYSTEM 






DESIDE POWER AND STEAM DEMANDS INPUT THE POWER 
AND STEAM DEMANDS INTO. THE CONTROLLER 




<?\ 3 


TRANSMIT THE POWER AND STEAM DEMANDS FROM 
THE CONTROLLER TO AN OPTIMIZING SYSTEM; 
MAKE AN OPTIMIZATION REQUEST FROM CONTROLLER 




< S 1 4 


SEARCH AN OPT 1 MAL SOLUTION 




1 5 


TRANSMIT THE OPTIM 
OPTIMIZING SYSTEM T 


AL SOLUTION FROM AN 1 
0 A CONTROLLER 1 




<S16 


CONTROLLER CONTROLS AN ENERGY 
SYSTEM OF THE PLANT 



FIG. 2 



17 



BNSDOCID: <EP 0731397A1TI_> 



EP 0 731 397 A1 




BNSDOCID: <EP 0731397A1TI_> 



EPO 731 397 A1 



# 



ELECTRIC POWER/ 
FLUID DEMAND 
INPUT SECTION 



2 5 



7 



2 1 



1 



FLNID OEMAND 
SYSTEM OPTIMIZING 

SECT ton 



22 



1 



LIST CONVERSION 
SECTION 



23 



LIST HOLDING 
SECTION 



FLUID INPUT SYSTEM 
OPTIMIZING SECTION 



2 4 



r 



FLUID GENERATION 
SYSTEM OPTIMIZING 
SECTION 



27 



r 



PLANT DATA 
HOLDING SECTION 



2 e 



LIST RETRIEVING 
SECTION 



2 8 



r 



OPTIMAL SOLUTION 
OUTPUT SECTION 



Z 9 



r 



FIG. 4 



BNSDOCID: <EP 0731397A1TI_> 



19 



EP0 731 397 A1 



( START ^ 




1 



S2 5 



S23^ 


YES 

f 




TCO][ j 0 0 0 0 | 


> 






f 


* — 


j = j -H | 


k = 0 | 



r S2 6 




✓S2 9 



TC0][j+k]-pi(j)+ P 2Ck) 
T[1][j+k][0]=J 

TCUCj+kHU-k 



S32 



k = k -h 1 



FIG, 5 



20 



BNSDOCID: <EP 0731397A1TI > 




T'tOJCi +k]-TCO]Cil+i>3(k) 

r [i ][j +n][o]-TCi)[j][o] 

TiUli +k]C1]»TCn[J][l) 
r [1Ki+k][«]-k 



k=k-H U 



FIG. 6 

21 

BNSDOCID: <EP 0731397A1TI_> 



EP0 731 397 A1 



■ b 1 1 AM AMUUM 

OF L3 


0 


r « ■ 


M1+M2+M4 | 


OPTIMAL STEAM 

DISTRIBUTION 

(SLS2,S3) 


(0, 0, 0) 


* «• p 


(Ml. M2, M4> ] 


OUTPUT POWER 


0 




P 1 (Ml) +p2 (M2) + p4 <M4) 



FIG. 7 



STEAM AMOUNT Of L4 


0 


Iff 


M3+M5 


OP 1 1 MAL STEAM 
DISTRIBUTION 
U3, S5) 


(0, 0) 


m ■ • 


0M3. M5) 


OUTPUT POWER 


0 


• l • 


p3 <M3) + p5 (M5) 



FIG. 8 



D 


0 


d 4 r 


uM 


STEAM AMOUNT 
OF L3 


0 


«l ► < 


Ml -S-M2+M4 


OPTIMAL STEAM 
DISTRIBUTION 


CO, 0. 0) 


« fc 4 


(Ml. M2, M4) 


OUTPUT POWER 


0 


► -1 ► 


pi (Ml) + P2 (M2) +P4 (MO 



FIG. 9 





0 


h J * 


vM 


•STEAM AMOUNT OF L4 


0 




M3+M5 


( S3, $5) 


(0, 0) 


* r * 


(M3, M5) 


OUTPUT POWER 


0 




p3 <M3) + p 5 <M5) 



FIG. 10 



22 

BNSDOCID: <EP 0731397A1TI_> 



* 



EPO 731 397 A1 



# 



u 


• Urn 




UM 


STEAM AMOUNT 

OF Li . I 


L 3 i 




L3 u 




CsmU &m2 ( sm3) 


1 4 t 


(sMI, sM2, sM4) 


OUTPUT POWER 


Em3 


h • a 


EM 3 



FIG. H 



V 


Vm 


■ ■ • 


VM 


STEAM AMOUNT OF L4 


L4 I 


^ ^ ■ 


L4 u 


OPTIMAL STEAM DISTRIBUTION 
( si, 55 J 


Csm3, sm5> 


r 1 1 . 


<sM3, sM5) 


OUTPUT POWER 


Em4 




EM. 4- 



FIG. 12 



TOTAL INPUT STEAM AMOUNT 
UN TERMS OF HIGH PRESSURE ) 


Um+Vin i 




UM-I-VM 


(STEAM OF L3, STEAM OF L4 ) 


(L3 Is L4f) 




(L3tu Liu) 


OPTIMAL STEAM DISTRIBUTION 
(SI, ,55) 


sm5) 


• • ■ 


(sM 1 

sMS) 


TOTAL OUTPUT POWER 


Em3+Em4 • 




EM3+EM4 



FIG. 14 



BNSDOCID: <EP 0731397A1TI_> 



23 



EP0 731 397 A1 



S5 1 



j =0 



S5 



j <J 1 + 1 2- [ ? 



1 



S55 



S53^ 


YES 


J-0 


j t[3][j]=-1 0 0 0 0 




< 










k-0 



'S5 & 




S59 



TC3JC j +k]«A[3][J]+BC3]tkl 
TC03[j -hk]=Aro][j]+SCO][k3 
TCULj +k3C0]^AClKJ] 
TC11LJ +k]Cl] = B[1][k] 
TC21C j + k]-AC2Kj]+BE2 3Ck] 



S6 I 



k = k + 1 



FIG. i 3 



24 



BNSDOCID: <EP 0731397A1TI_> 




EP 0 731 397 A1 





Um+Vm 




UM+VM 


(STEAM OF L3, STEAM OF L4 ] 


(L3K Hi) 




CL3u, L4u) 


STEAM DISTRIBUTION 

(S|,.\. ( S5) 


(sml 

sm5) 




CsMI 

sMS) 


TOTAL OUTPUT POWER 


Em3+Em4 


i 4 J 


EM34-EM4 


FUEL DISTRIBUTION 


(*, *, *, *) 




(*, *, *, *) 


TOTAL FUEL COST 






* 



FIG. 15 



TOTAL INPUT STEAM AMOUNT , 
( IN TERMS OF HIGH PRESSURE! 


Um+Vm 




UM+VM 


( STEAM OF L3 , STEAM OF L4 ) 


(L3 fl, L4 i) 




CL3 Us L4 u) 


STEAM DISTRIBUTION 
(SI,... ,s5 \ 


(sml , . . . , 

SfTlS) 




CsM 1 , , . . , 

sM5) 


TOTAL OUTPUT POWER 


Em3 + Em4 


4 4 4 


EM 3 +EM4 


FUEL DISTRIBUTION 






C*, *, *, *) 


TOTAL FUEL COST 








TOTAL COST 




i ► ► 


* 



FIG. 16 



25 

BNSDOCID: <EP 0731397A1TI_> 



EP0 731 397 A1 



<2 

o 55 



* - i 

J. S5i-32 
po- OSuj 



T — 



Sill 



CP 



CM 



cr 

o 



CM ^ 



O 



CD 



t>— -JO 





26 



BNSDOCID: <EP 0731397A1TL> 



# 



EP0 731 397 A1 



* 




BNSDOCID: <EP 0731397A1TI_> 



27 



EP0 731 397 A1 




FIG. 19 



BNSDOCID: <EP 0731397A1TI_> 



2B 




F I G. 20 FIG. Z\ 



FIG. 22 



FIG. 23 




BNSDOCID: <EP 0731397A1TI_> 



29 



EP0 731 397 A1 




GLOBAL OPTIMAL POINT 
(TRUE OPTIMAL POINT) 



F I 6, 24 



30 

BNSDOCID: <EP 0731397A1TI_> 



♦ 



EP0 731 397 A1 



* 




FIG. 25 



PLANT MODEL 

CHARACTERISTIC 

TURBIN OPERATION 
RESTRICTIONS 

BOILER 

CHARACTERISTIC 

BOILER OPERATION 
RESTRICTIONS 
, STEAM BALANCE efc 



PLANT STATE 
ELECTRIC LOAD AMOUNT 
STEAM LOAD AMOUNT 
ENERGY BUYING UNIT PRICE 



31 
<± 

TOTAL ENERGY 
COST 

MINIMIZING POINT 

CALCULATION 
SECTION 



[ TUNING' 1 

i REQUIRED I 

I CALCULATION j 
L PARAMETERS 



FIG, 26 



OPTIMAL . 
SOLUTION 



BNSDOCID: <EP .0731397A1TI_> 



31 



* 



THIS PAGE BUkUK (uspto) 



